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*>— ^ ' Abstract. We present a class of models in which the primordial metric fluctuations do not necessarily obey Gaussian statistics. 

^^ ' These models are realizations of mechanisms in which non-Gaussianity is first generated by a light scalar field and then 

O transferred into curvature fluctuations during or at the end of inflation. For this class of models we present generic results for 

^ ■ the probability distribution functions of the metric perturbation at the end of inflation. It is stressed that finite volume effects 

^O ' can induce non trivial effects that we sketch. 
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> 1. INTRODUCTION 

O. 

C^ ■ Cosmic Microwave Background (CMB) observations offer a precious window for the physics of the early Universe. 

And it is well known that inflation generically predicts Gaussian initial metric fluctuations with an almost scale 
invariant power spectrum [1]. With the advent of large scale structure and CMB surveys however it will be possible 
\^^ ' to test in details the statistical properties of the initial conditions. So far no non-Gaussian signal has been detected in 

f^ . CMB data [2] but the number of modes that can be probed is still small. In large-scale structure surveys the number 

of independent modes that are observed is large but the difficulty arises from the non-linear gravitational dynamics [3] 
which can shadow the primordial Non-Gaussianity (NG). Having at our disposal models of inflation in which non- 
Gaussian adiabatic metric fluctuations are generated can then serve as a guideline for designing strategies for detecting 
O \ primordial non-Gaussianities, in particular in the context of the Planck mission. 

4— > ' A number of propositions have been made which circumvent the usual bounds on primordial NG generic inflation 

^ \ (e.g. slow-roll single field inflation) predicts (see below). This is the case in particular when primordial adiabatic fluc- 

tuations are subdominant compared to density fluctuations generated during the pre- or re-heating phases. Examples 
of such models are the curvaton model (which relies on assumptions on the reheating phase [5]) or recent propositions 
of multiple-field inflation which are based on relatively well understood phases of preheating [6] . 



2. GENERIC INFLATION 

As mentioned in the introduction, we are interested in models that can produce sufficiently large non-Gaussianity. 
To be more precise, that should be the case for at least the range of modes that are observationally relevant, i.e. that 
corresponds to the large scale structure scales. 

Let us start by considering a single field, (p, in slow-roll inflation. The Klein-Gordon equation for its perturbation, 
5(p, is of the form 

5(? + 3H5(?-^5(? = -V"5(?-V"'-^ + ... (1) 

a^ 2 

During the slow-roll regime, |//| ^ H^ and (p <c 3// 9 so that' 

3i/^M| ~ y, //^ ~ -V' (2) 



We tiave set M.'- = ?,7tG. 
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FIGURE 1. The trajectories of the fields in the plane ((p-X), once smoothed on a scale R. Before horizon crossing (R < H), 
the trajectories behave quantumly because the quantum fluctuations are active up to scale H and are not smoothed out. After the 
horizon crossing (R > H), the trajectories can be treated as classical trajectories. Note that this transition happens at different time 
for different values of R. The bundle of classical trajectories then evolves in the two dimensional potential and its cross section 
evolves with time. In the case sketched here, because of the bending of the potential valley, the isocurvature modes induce metric 
fluctuations because of differences in length of the trajectories. The net result is a transfer of modes between isocurvature and 
adiabatic directions. 



and the slow-roll conditions can be expressed as e <C 1 and | Tj | <C 1 with 

V' .V" 



(3) 



In order for the fluctuations of the scalar field to be large enough, one needs the mass of the field to be much smaller 
than H, in which case one gets {5(p) ^ H. Furthermore to get the correct amplitude for the primordial fluctuations, 
one should have V^I^/{V'mI) - 10"^ so that 



//-10"^eM4. 



(4) 



The range of wavelengths that correspond to the observed large scale structure exits the Hubble radius over the number 
of e-foldings 

Nx=HM. (5) 

During that period, the slow-roll parameter 7] has varied over the range At] ^ {t, — rie)A(p/M4, with ^ = M^iV'" /V). 
From Eq. (4), one deduces that 

A(P/Ma eHAt. (6) 

Now, if the quadratic term in the rh.s. of Eq (1) is dominant over the linear term then V" /V'" < 5(p ^ H and, using 
Eq (4), one deduces that Arj > lO^TjA^;^. This will induce a rapid breakdown of the slow-roll inflation. This can be 
understood from the fact that the potential has to be both flat enough for the fluctuations to develop and steep enough 
for the non-linear terms not to be negligible^ 



3. MULTIPLE FIELD INFLATION 



The previous analysis shows that it is not possible to get significant NG effects over the observable wavelengths with 
single field inflation. To obtain such effects, extra fields are necessary. A class of multiple field inflationary model 
corresponds to a Lagrangian where A^ scalar fields are minimally coupled to the metric, e.g.. 
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Y^d^(Pjd^(pj-V{(pi,...,(pN), 



(7) 



■ j=i 



^ A way round to this argument is to consider potential with a sharp feature [4] but in that case the non-Gaussianity is associated with a departure 
from scale invariance and is located on a very small band of wavelengths. 



where R is the Ricci scalar. It follows that the background Einstein equations take the form, 

//2^^-l ^(p2+2yj, //=-47rG^(p2, (pj + 3H,Pj = -Vj, (8) 

where a dot refers to a derivation with respect to the cosmic time, t, H = cija and V, refers to a derivation with respect 

to (?j. 

Then, let us consider the simplest case in which we have only two scalar fields ^i and ^2- We decompose these two 
fields into the genuine inflationary direction 5^ and a transverse direction %, taking advantage of the direction defined 
by the inflationary direction. 



X ) \92 

the angle Q being defined by 

cose = ^ , sine = ^ . (10) 

From these definitions it follows that % is stationary along the classical inflationary trajectory. To be more precise, the 
Klein-Gordon equations imply that. 
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8(p + 3H8(p-^ = ---^8(p-^^^x, (11) 



This set of equation clearly opens the way to new phenomena. Indeed, nothing prevents the fluctuations in the x 
directions to develop significant non-Gaussian properties. This would be the case for instance for generic quartic 
potentials. Those field fluctuations can appear as metric fluctuations if, at some stage before the end of the inflationary 
period, the x fluctuations can be transferred into the adiabatic fluctuation. According to Eq. (11), this is possible 
provided d^V / {d 8<f>d x) is not zero. Geometrically speaking, it corresponds to a bent in the inflationary trajectory^ [7, 
8]. This idea is sketched on diagram 1. 



3.1. A simple working example: extended hybrid inflation 

Although generic predictions can be made for such a class of models, it is to be stressed that such a scenario does not 
require elaborate or exotic field theory constructions. Let's for instance consider the following multiple-field potential 
presented in [10], 

V((/),Z,CJ) = im2</)2 + A;^4^|((j2-cJo2)' + |cj2((/)cosa + xsina)2. (13) 

It is an extension of the hybrid model with one extra field. Here, one field, 0, is the inflaton; the second field is a 
light scalar, j, with quartic coupling < A < 1 and the third field, a, is coupled to the two others so that the end of 
inflation is triggered when a undergoes a phase transition and cTq is the final vev of a. One extra parameter a describes 
the mixing angle between and x in their coupling to a. This parameter will determine the ratio between the initial 
isocurvature and adiabatic modes in the final metric fluctuations. Note that a quartic interaction has been introduced 
since this is the only one that does not require further fine tuned parameter (any other polynomial type interaction 
would require dimensioned parameter tuned to H, see [8]). 

In this model, as for hybrid inflation, the inflationary stage ends when the effective mass of a vanishes, that is when 

g(0cosa + xsina)^-2jU(7o ==0. (14) 



Note that this bent can also be encoutered during the preheating phase. An example of such a model in presented in [9]. 



The value of (j) at the end of inflation is therefore 



_ ±y/2n/gao-xsina 
COS a 

For ^ > ^end^ G — and the two fields evolve independently: (p drives the inflation while x develops non-Gaussianity. 
The amount of non-Gaussianity of x then depends only on X and on the total number of e-foldings between horizon 
crossing and the end of inflation. 



3.2. Test quantum scalar field in de Sitter space 

The non-Gaussian properties of the x field during the inflationary phase are those developed by a self interacting 
scalar field in an expanding universe. The full resolution of this problem in quite involved. It can be addressed with a 
Perturbation Theory approach that can be done either at a quantum level or, once the horizon has been crossed, at a 
classical level (assuming the scalar field behaves like a classical stochastic field). 

In [12] we present the result of the computation of the high order correlation function of a test scalar field in a 
de Sitter or quasi de Sitter background. The calculation is based on an expansion of the field evolution from the free 
field solution as summarized in the following. For a minimally coupled free quantum field, the solution can indeed be 
decomposed in plane waves as. 



vo(x,77)= d\ vo(fc,77)^ke"''' + vS(^,J7)^: 



/k-x I ,,*/'/, «\/^'''a~'''*^''^ 



(16) 



where we have introduced v = a 2, a hat referring to a quantum operator and fe^ and fo^ and the annihilation and creation 
operators for a particle of momentum k in a Bunch-Davies vacuum. In the massless limit, the free field solution is 



1 \ e-i*^") 



One can then express perturbatively the A^-point correlation functions of the interacting field, j, in terms of those of 
the free scalar field. The equal time correlators are expectation values of product of field operators for the current time 
vacuum state. Such computations can be performed following general principles of quantum field calculations [11]. 
The simplest formulation is to apply the evolution operator U{rjo,rj) backward in time to transform the interacting 
field vacuum into the free field vacuum at an arbitrarily early time tjq so that, 

(vu,...vu„) = (0|f/-'(77o,T7)vk,...vk„f/(77o,77)|0) (18) 

where |0) is here the free field vacuum "*. 

The evolution operator U can be written in terms of the interaction Hamiltonian, Hj, as 

f/('?o,'?) = ^exp('-i|''d77'///(77')') (19) 

where ^ is the time ordered product operator^. Eventually the connected part of the above ensemble average at a time 
rj reads, 

(vk. . . . vk„)c = -i Pdri' (0| [vk, • • • vk,.,H/(77')] |0), (20) 



where the brackets stand for the commutator. 



** It is to be noted that these calculations do not correspond to those of diffusion amplitudes of some interaction processes in a de Sitter space 
(see [13] for a comprehensive presentation of those calculations.). When one tries to do these latter calculations with a path integral formulations, 
mathematical divergences are encountered as it has been stressed in [14, 8]. With that respect de Sitter space strongly differs from Minkowski 
space- time. 

^ Then the inverse operator for U reads, [/"'(jjo.'j) = ^exp (+ijn dJ]' //;(j]') 1 where .T is the inverse .T product. 
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FIGURE 2. Behaviour of the function Q4 as of function of time. The transition to the superhorizon behavior (dashed line) is 
shown. The function Q4 is shown here for a "square" configuration (fci = k2 = kj = ^4) as a function of ktrj = ^fcjT]. 

The result for the four-point function in case of a quartic interaction (with parameter X) can finally be written as, 



H^ 






where the vertex value, V3, reads 



^3(M) = ^^[C('?,fcO+iogHEfcO] 



(21) 



(22) 



In this expression ^{rj,ki) can be expressed^ in terms of the four wavelength A:,. It is finite for any values of 77 and k, 
(in particular in the super-Hubble limit). 

When the term log (— Tj L^O is large (and negative), that is when the number of e-foldings, Ne, between the time of 
horizon crossing for the modes we are interested in and the end of inflation is large, the vertex value is simply given 

by, 

V3(M) = -AA^,/(3//2), 

which corresponds exactly to what a classical stochastic approach would give [8]. On Fig. 2 one can then appreciate 
the transition from a regime which is dominated by quantum fluctuations to a regime following a classical evolution. 
In this figure, what is plotted is the reduced four-point correlation function, Q4, defined as. 



G4({fc,}) 



P4{kxMMM) 



P2{ki)P2{k2)P2{k3)+sym. 



(23) 



where P2{k) is the free field power spectrum. 

These results illustrate one quantitative aspects of the NG properties that self-coupled scalar fields can develop, 
namely high order correlation functions, whose computation is given here at tree order. These calculations are on solid 
ground. They do not give however a complete prescription for the description of the statistical properties of the field. 

3.3. PDF and finite volume effects, a classical stochastic approach 

The use of a classical perturbative approach can provide us with a complete description of the expected properties 
of the resulting field properties. In particular it is possible to infer the whole one point probability distribution function 
(PDF) of the field value. 



note that f differs from the one introduced in [12]. 



The approach detailed in [8] and [15] is based on the expansion the filtered field in terms of the coupling constant 



as 



_v(0) 



Zs('7)-Zr('7) + Z'n'7) + 



1)/ 



(24) 



The subscript S stands for the fact that the field has been convolved with a given smoothing window function in such 
a way that only large enough scales are taken into account. Then y '*'' represents the value of the filtered field when 
the self-interacting term in the potential is dropped. In the slow-roll regime the field x^ follows the Klein-Gordon 
equation (12). Then xf is simply constant during the slow roll evolution. It represents the initial conditions for the 
non-linear field evolution. At first order term in X, Xs = Xs ' evolves according to 



3i/is"'- 



X_ 
3! 



(0) 



(25) 



for a quartic potential. In this approach the treatment of the filtering of the r.h.s. of this equation is very crude. This 
is the price to pay for using this approach. The results we are going to find can anyway be checked against the more 
rigorous calculations based on the computation shape of the tri-spectrum presented in the previous subsection. 

Resolving Eq. (25), leads to the expression of jj'^' as a function of j'"'. Xs is actually built out from modes whose 
wavelengths are larger than the smoothing scale. During inflation, Xs can be viewed as a classical random field as 
soon as the smoothing scale has crossed the comoving horizon. Today, Xs is made of a the superposition of modes 
that are super-horizon and of modes which are now within our horizon. It is important to understand that, although 
the ensemble average of Xs actually vanishes, its geometrical mean at the survey size, as given by observations, is 
expected to be non-zero (since at least super-horizon modes are out of reach). 

The value of Xs at our horizon scale, say Xs > appears therefore as a new free parameter. Its value is determined by 
the peculiar realization of the inflationary path our Universe has followed. In our model case xi is bounded (A > 0) 
and its range can be estimated through a Fokker-Planck equation [16, 15]. Phenomenologically it means that one 
should decompose r '"' as. 



,(0) 



Sx, 
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(26) 



where dx^^ is then Gaussian distributed with zero average and given variance (Jg. Similarly, 



Xs = SXs+X 



(27) 



where x corresponds to the averaged value of x at our horizon size. 
The equation of evolution (25) can be solved to get 



z(')(f) = -A(f-f„) 



X. 



(0) 



18// 



(28) 



which also reads 



Xi'\t) 



XN, 



18//2 V^s 



(0) 



Ne being the number of e-foldings between tn and the end of inflation. 
These results imply that 



X ^ Xi 



(0) 



XNe 

18//2 



(0) 



-^x 



(0)^2 



(29) 



(30) 



which explicitly shows that x and x& are equal at leading order in A. It also gives 
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SXs ^ ^x 



'"'-l^lNI^-3 



(5zf) -. 



x+^sxi'^f 



(31) 



It is straightforward to see that dx^ has acquired a non zero third order moment. 
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(32) 
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FIGURE 3. PDF of 3%^ = hx^ — % for different values of %. The dashed line corresponds to a Gaussian distribution; the dot- 
dashed line to the deformed distribution of 8%^ when XNe/tfi = 1 and ^ = and the solid lines when the latter equals 0.5 and 
1. 

at leading order in X. This is a finite volume effect in the sense that it exists for a fixed (not ensemble averaged) value 
of J. This effect cannot a priori be neglected. It can be shown that X should be of the order of H /X ''^, which implies 
that the reduced skewness of Sjj, ^xlK^xTf' ^ X^'^NeOg/H is significant as soon as X^'^^Ng approaches unity. 
An attentive reader cannot but notice that the evolution equation for j„ can in fact be solved in 



X, 



z, 



(0) 



(33) 



9H^ 



Xi 



(0) 



and the distribution of Xs can then be infered from that of x^ ' assuming the latter is Gaussian distributed with a 
non-zero mean value '' . 

In figure 3, we present the deformation of the PDF of 5x^ while x is varied. As expected, it shows that when x 
is not zero, the PDF gets skewed in a way that can be easily understood: when x is positive it gets more difficult to 
have excursion towards larger value of j^, but easier to roll down to smaller values. It is as if the field x was actually 
evolving in the potential A(x + x)/4!. 

What should then be expected for observations? In this family of models, the surviving couplings in the metric, 
are expected to be, to a good approximation, equivalent to those induced by the superposition of two stochastically 
independent fields, a Gaussian one and one obtained by a non-linear transform of a Gaussian field with the same 
spectrum, e.g. Eq. (33). In other words the local (Bardeen) potential would read. 



<I>(x) =:cosa 4>i(x) +sina ^[4>2(x)] 



(34) 



where ratio between the initial adiabatic and isocurvature fluctuations is described by a mixing angle a and where the 
function ^ takes into account the self coupling of the field x that gave rise to this part of the metric fluctuations. This 
description is the starting point of investigations for the CMB properties presented in [17]. 

4. CONCLUSIONS 

We have shown here that it is possible to build explicit models of inflation where primordial metric fluctuations can 
develop significant non-Gaussian features. The class of models we have identified induce NG properties that can 



As noted in Ref. [8], such a simple variable change implicitly incorporates "loop order" effects that, because of sub-Hubble physics, are not 
necessarily correctly estimated. In that paper we developed a more elaborated method which allows the reconstruction of the PDF from the only 
tree order contributions of each cumulant. As the two approaches eventually give the same qualitative results we restrict here our analysis to the 
simplest method. 



be described by a few identified parameters. Two of them are of fundamental origin, the amphtude of the quartic 
coupling and the mixing angle; a third is due to finite volume effects and its value depends on the peculiar realization 
of the universe we live in, the transverse field value at Hubble size. Unlike most models proposed in the literature, 
they provide us with well defined properties of the metric fluctuations that can serve as a starting point for further 
investigations on its various observational aspects, both for CMB and large scale structure surveys [17]. In particular, 
the models do not require to deal with the more intricate second order dynamics. 
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